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Abstract 

We present a model for the kinetics of spontaneous membrane domain (raft) assembly that 
includes the effect of membrane recycling ubiquitous in living cells. We show that the domains 
have a broad power-law distribution with an average radius that scales with the 1/4 power of the 
domain lifetime when the line tension at the domain edges is large. For biologically reasonable 
recycling and diffusion rates the average domain radius is in the tens of nm range, consistent 
with observations. This represents one possible link between signaling (involving rafts) and traffic 
(recycling) in cells. Finally, we present evidence that suggests that the average raft size may be 
the same for all scale-free recycling schemes. 
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Liquid-liquid phase separation in mixed membranes is now a well known phenomenon. 
Separated domains in vesicles have been observed by fluorescence labelling and the size 
of circular sphingomyelin (SM)-enriched domains can reach almost half that of a 20/im- 
sized vesicle These domains quickly reassumed their circular shape if deformed, 

and merged with each other to create circular domains, all phenomena consistent with the 
existence of strongly immiscible liquid domains. There has been much recent biological 
interest in "lipid rafts" which are believed to be liquid ordered membrane microdomains 
containing certain proteins and enriched in glycosphingolipids and cholesterol DDQ . These 
are thought to have many important functions, including in signal transduction and in the 
sorting of proteins. Experimental estimates of the size of lipid rafts in vivo are in the few 



tens to 100 nm range 
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121. We refer to these as "intermediate" 



or even smaller 

sized domains in what follows to distinguish them from the micron scale (or larger) domains 
observed in reconstituted membranes that are close to equilibrium. 

We first review why intermediate sized domains are not expected to form on two cont- 
inent membranes at equilibrium in the absence of long range interactions j3] or recycling 
l^j ]. We assume that the surface coverage of domains remains small and write a simple 
Flory-Huggins model for the grand potential energy per area as a function of the distri- 
bution of the sizes of domains, defined by the dimensionless concentration c n of domains 
containing n monomers. All concentrations are per monomeric area s = ir(b/2) 2 with the 
monomer diameter perhaps b ~ 5nm if we identify the effective monomeric unit in vivo as a 
typical raft-resident proteins together with its associated lipid "skirt" [lfi| . 

oo 

^ = Cn ( log ~ ~ ^ n + 7v/ ™) (!) 

n=l 

Here and below all energies are measured in units of k B T. A line tension 7/(2 \j it s ) acts at 
the perimeters of all the domains. The chemical potential \x is conjugate to the total area 
fraction of domains <fi = X^Li n c n- Thermodynamic equilibrium then corresponds to 
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The average domain size is n = <j)/N with iV = XmLi c « ^ e total density of domains of all 
sizes. When is very small ~ log0 is large and negative and almost all domains are 
very small (monomeric) in size. For \i = there exists a critical total area fraction <p c = 
Ylri=i ^e~ 7v/ ™, again made up of small domains provided 7 3> 1. There is a phase transition 



for > C beyond which the area fraction of finite-sized, essentially monomeric, domains 
<ft c remains constant but the distribution becomes highly bimodal with a few very large 
domains accounting for the remainder <fi — C [22]. This is at odds with many experimental 
observations and hence motivates a non-equilibrium model, involving recycling, which will 
occupy the remainder of this paper. 



NON-EQUILIBRIUM: RECYCLING 

The evolution of the domain size distribution on an infinite membrane may be written 
via the following master equation 



dc n 
~dt 

1 ,pn-l if 7 /o\ 

~ 2 ' J Tfi — 1 m n—m { ^n—m L "m n m,n—m u n \ u j 



Here a{n) controls the lipid recycling, as will be discussed in further detail below; o = 
when there is no recycling. The kernals k n , m and k' n m control, respectively, the rates of 
domain scission in which one domain of n + m monomers breaks into two, of size n and 
m, and domain fusion, in which two domains containing n and m monomers fuse to form a 
single domain of size n + m, see Fig la. A similar approach has been rather successful in 
describing the kinetics of wormlike micelles 

We assume that two domains fuse whenever they diffuse into contact. Thus k! = D/s 
with D a characteristic diffusion constant for the domains and s the area of the smallest 
(monomeric) domains. The characteristic, microscopic "diffusion" time scale t d = 1/k' = 
s/D m 10 _5 s for biological membranes. We propose to neglect an y si ze-dependence of the 
diffusion constant of the domains, which is a fair approximation [l7(. This simplifies the 
analytic analysis and should only weakly affect our results, which rely primarily on the fact 
that the recycling rate is much slower than l/r D . 

The scission kernal k n ^ m is now determined by the principle of detailed balance, an ap- 
proach that is appropriate provided the longest intra-domain relaxation time following a 
fusion or scission event is shorter than the domain collision time. This should be satisfied 
on average for fluid domains where the rate of subsequent inter-domain events should 
then converge to that found at equilibrium for each (pair of) domain(s), in spite of the fact 
that the size distribution may be far from equilibrium. By inspection of Eq (jHJ) together 
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with Eq © 

where k' nm = k' = l/r D is treated as a constant in what follows. This kinetic scheme, 
without recycling^ will yield asymptotic domain growth reminiscent of the coarsening of 
crystal domains |l9j ]. 

In order to analyse the effect of lipid recycling we first consider the "monomer deposition 
/ raft removal" (MDRR) recycling scheme in which raft lipids and proteins are brought to 
the membrane at random as single 'monomer' sized units with a rate j on and entire rafts 
are lost from the membrane at random with a rate j a g, irrespective of their size, see Fig lb. 
Thus <j(n) appearing in Eq (JHJ) takes the form 

= jonS n ,l ~ joSC n (5) 

with Sij the usual Kronecker delta. It is easily shown that 6 — jon/jos at steady state. 
In general Eq (JHJ) requires numerical solution, see Fig 2 and |23j ]. although an asymptotic 
solution is possible in the most interesting regime 7 3> 1, as discussed below. 

It can be seen from Fig 2 that the domain size distribution is broad, indeed there is 
significant contribution to the total area fraction from domains with n < n 2 . The distribution 
depends only weakly on 7 when 7 is itself large because all domain scission events are then 
rare. 

To investigate this 'scissionless' large 7 regime further we proceed by neglecting all scission 
terms in Eq Q and obtain 

dc 1 

_ ^ = a (n) - k'N c n + k' - c n _ m c m (6) 

where we will later have to check our solutions for self-consistency, which will translate to 
establishing a lower bound for 7. 

The equations c n = for n — 1, 2 . . . can then be used to build up c n recursively. For 
MDRR recycling given by Eq Q this yields 

C n = A n (jonr D )7(joffT D + iV) 2 ™" 1 (7) 

with A n = (2n — 2)!/(2 n_1 n!(?i — 1)!). Eq ((Zj) still involves iV and so the size distribution 
is only explicitly determined after solving for iV = c n . The resulting c n appears as the 
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solid curve on Fig 2. In the most biologically relevant regime |2J] 4>/(j aT D ) ^> 1 and 
N = ^2(fij fiT D . Eq[7|then approaches 



6j off r D /(2vr) n^' 2 for 1 « n < n 2 (8) 
which appears as the dashed curve on Fig 2. The average domain area in this regime is 



n = v/0/(2 Joff r D ) (9) 

as shown in Fig 3 (dashed line). 

When scission is rare the average domain radius is found to be R = 10-70nm, for 6=1- 
5nm respectively for recycling with rate j Q s = 10~ 2 s _1 , see Fig 3. As can be seen from 
the raft area Eq © the raft radius scales only with the 1/4 power of the recycling rate or 
diffusion time and hence is in the tens to ~ lOOnm range for all reasonable physiological 
values. Restoring dimensionality the threshold line tension 7* « 10 is about 0.6/c B T/nm 
when 6 = 5nm. The line tension between different membrane phases may have several 
origins. Hydrophobic thickness mismatch [25] gives rise to line tensions that depend upon 
the elastic properties of the membrane and are in the fc B T/nm range. Other, more directly 
chemical or ordering-related, incompatibilities may also contribute. Thus the asymptotic 
results of this section, in which we assumed that the domains rarely break, may have wide 
applicability and full numerical solution of the master equation may not always be necessary. 

The scission rate Eq has a maximum at m = 1 corresponding to shedding single 
monomers. For domains with n» 1, such as those of the average size, the maximal rate 
is A;i, n -i ~ e'^k' which can be used to establish a self-consistency condition: Only if the 
number of monomeric scission events in the typical residence time of a domain is much less 
than n will the scission rate be negligible. The product of this rate and the lifetime j~Q is 
much less than n whenever 7 3> 7* with 7* = \ log -r^- — . 

Another intuitive 'scale free' recycling scheme is the "monomer deposition / monomer 
removal" (MDMR) scheme in which raft lipids and proteins are again assumed to be brought 
to the membrane at random as small 'monomer' sized units with a rate j on but domain 
material (raft lipids and proteins) are removed from the membrane in monomeric units, 
at random, with a rate j ff , irrespective of the size of the raft on which they reside. While 
perhaps of less biological relevance this scheme, together with MDRR, form the most extreme 
examples of an entire class of possible recycling schemes: they correspond to the removal 



of monomers, through to whole rafts, respectively and involve no characteristic size scales 
for the recycling. As such it would be surprising if any intermediate 'partial raft' recycling 
scheme produced behavior that was not bounded qualitatively by the two extremes that we 
consider. Within this scheme 



with = jon/joff at steady state, as before. The general solution of Eq (JHJ) is again obtained 
numerically, see Fig 4. 

Fig 5 shows that the mean domain size is the same for both of these recycling schemes, 
suggesting that this is true for all scale-free recycling schemes. To understand this note 
that the lifetime of every domain monomer is the same under both schemes. It is possible to 
propose other schemes that are not scale-free but rather have intrinsic size scales and it is far 
from clear that such schemes may not be biologically relevant. Such scales might correspond 
to a characteristic raft size at which endocytosis occurs. While it is straightforward to treat 
such schemes by an appropriate choice of crin) the mechanisms by which lipids are recycled in 
the cell are still rather poorly understood [20] and a thorough investigation of the numerous 
possible schemes lies outside the scope of the present work. 

DISCUSSION 

We have shown how the non-equilibrium nature of membranes can lead to steady-state 
domain sizes that are intermediate in size, typically in the tens to lOOnm range for all bio- 
logically reasonable recycling rates and membrane diffusion constants. This result seems to 
agree well with experimental observations of lipid rafts on the plasma membranes of different 
cells and is marked contrast to the large domains observed on artificial membranes as they 
approach equilibrium. We believe that our results will be of interest to those working to re- 
alise membrane recycling in model (bio) chemical systems as well as those seeking candidates 
models that may provide a better understanding of lipid rafts in living cells. 

It is now being realized that signaling and traffic in cells may be closely related processes. 
Our model shows how the regulation of membrane traffic (recycling) might simultaneously 
control the raft sizes. It is also quite plausible that the size of rafts controls certain aspects 
of their function and this would give a direct connection between signaling and traffic. 




(10) 
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We are currently studying the potentially important effects of finite (cell) membrane 
area, which may result in a failure of the mean field approach for the largest rafts. Also of 
interest is the response of the raft size distribution to perturbations in the recycling, such 
as a step change in the monomer deposition rate. This will provide a biologically accessible 
way of relating response time(s) to this model and hence the rate of the underlying recycling 
mechanism. 
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FIGURE CAPTIONS 



Fig 1. The model: (a) all domains can undergo scission and fusion (b) domains are 
deposited and removed from the membrane stochastically. 

Fig 2. The steady state domain size distribution for MDRR recycling with = jon/jos = 0.1, 
JoffT" D = 10~ 3 and 7 = 2 (o) 4 (A) and 6 (□). The mean domain contains 1.3, 2.3 and 5.3 
monomers respectively, rather small because a fast recycling rate was chosen for numerical 
convenience. The asymptotic variation for large line tensions Eq (|ZJ) is shown (solid line) as 
is Eq (JSJ) (dashed line) appropriate for sizes 1 << n « n 2 . 

Fig 3. The variation of the mean domain size n and radius R with <fi/(j gT D ) (solid 
line) assuming that the domain scission rate is negligible and b = 5nm. The power law 
variation for large, but physiologically relevant, (p/(j sr D ) is also shown (dashed line) as is 
the average residence time of the domains t s, with r D = 10~ 5 s and = 0.1. 

Fig 4. The steady state domain size distribution for the MDMR recycling scheme. 
The same values as Fig 2 are used except 7 = 3 (o) 6 (A) and 9 (□). The mean domain 
contains 1.5, 5.2 and 7.0 monomers respectively, larger than in Fig 2 because of the larger 
line tension 7. The asymptotic behavior for the MDRR scheme at large line tensions Eq (J2J) 
is shown (solid line), as is Eq (jHl (dashed line), but it overestimates the very small number 
density of the largest domains present under this scheme. 

Fig 5. The steady state mean domain size n as a function of 7 for both MDRR 
(polygons) and MDMR (crosses) recycling, with = 0.1 and j s T ly — 10 _1 (x and A) 
joffT D = 10~ 3 (+ and □) j ffT D = 10~ 4 (crosses and pentagons). The asymptotic values of 
n for 7 ^> 7* from Fig 3 are shown as dashed lines at n — 1.4 (j ffT D = 10 _1 ), n = 7.6 
(joffT D = 10~ 3 ) and n = 23 (j sT D = 10~ 4 ). The corresponding values of the crossover 
tension 7* are 2.6, 5.0 and 6.1 respectively and are close to each inflection point. 
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